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ABSTRACT 

It seems to be feasible in the near future to exploit the properties of left-handed metamaterials in the telecom or 
even in the optical regime. Recently, split ring-resonators (SRR's) have been realized experimentally in the near 
infrared (NIR) and optical regime.^' ^ In this contribution we numerically investigate light propagation through 
an array of metallic SRR's in the NIR and optical regime and compare our results to experimental results. 

We find numerical solutions to the time-harmonic Maxwell's equations by using advanced finite-element- 
methods (FEM). The geometry of the problem is discretized with unstructured tetrahedral meshes. Higher 
order, vectorial elements (edge elements) are used as ansatz functions. Transparent boundary conditions (a 
modified PML method'^) and periodic boundary conditions^ are implemented, which allow to treat light scattering 
problems off periodic structures. 

This simulation tool enables us to obtain transmission and reflection spectra of plane waves which are in- 
cident onto the SRR array under arbitrary angles of incidence, with arbitrary polarization, and with arbitrary 
wavelength-dependencies of the permittivity tensor. We compare the computed spectra to experimental results 
and investigate resonances of the system. 

Keywords: Left-handed metamaterials, split-ring resonators, finite-element method, Maxwell's equations, trans- 
parent boundary conditions, periodic boundary conditions 
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1. INTRODUCTION 



With the advances in nanostrncture physics it has become possible to manipulate light on a lengthscale smaller 
than optical wavelengths.^'^ It is now possible to construct periodic structures made of artifical nanostructures 
with macroscopic properties which do not occur in nature. These nanostructures are large on the atomic scale, 
therefore they can be of complex geometry. But they arc small on the scale of the wavelength of the illuminating 
light, therefore the system has properties of an effective homogeneous material, in particular a macroscopical 
electric permittivity, e, and magnetic permeability, /i. 

Of special interest are metamaterials with both, negative e and negative //.^ This leads to a negative index 
of refraction^' and allows in principle to overcome limits in the resolution of optical imaging systems.^ 



2. ARRAYS OF SPLIT-RING RESONATORS 

Split-ring resonators can be understood as small LC circuits consisting of an inductance L and a capacitance 
C . The circuit can be driven by applying external electromagnetic fields. Near the resonance frequency of the 
LC— oscillator the induced current can lead to a magnetic field opposing the external magnetic field. When the 
SRR's are small enough and closely packed - such that the system can be described as an effective medium - 
the induced opposing magnetic field corresponds to an effective negative permeability, < 0, of the medium. 

Arrays of SRR's with resonances in the NIR and in the optical regime have been fabricated using electron- 
beam lithography on a 1mm thick glass substrate coated with a 5nm thick film of indium-tin-oxide (ITO). 
Figure 1 shows electron micrographs of a produced sample. The 'U'-shaped SRR's are produced with a thickness 
of A2; = 30 nm. Details on the production can be found in previous works. ^'^ 
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Figure 1. Electron micrographs of an array of SRR's with different magnifications. Field of view corresponds to approx- 
imately 8.5/umx7.6/tim (a), resp. 1.25 /xmx 1.0 pm (b). The total size of the fabricated array is lOO/xmxlOO/im. (See 
original publication for images with higher resolution.) 

Due to the small dimensions of the LC circuits their resonances are in the NIR and optical regime.^ The 

elementary vectors of the periodic array, ai = (315 nm, 0, 0) and a2 = (0, 330 nm, 0), are smaller than NIR and 
optical wavelengths, therefore in these regimes the system is well described as an effective medium. 

In what follows we investigate numerically the transmission of light under oblique incidence through these 
samples. 



3. LIGHT PROPAGATION IN PERIODIC ARRAYS OF SRRS 

We consider light scattering off a system which is periodic in the x— and y— directions and is enclosed by 
homogeneous substrate (at Zsub) and superstrate (at Zgup) which are infinite in the — , resp. -|-^;— direction. Light 



Figure 2. Geometry of the unit cell in the array of SRRs. The surrounding air and the substrate are not indicated. (See 
original publication for images with higher resolution.) 



propagation in this system is governed by Maxwell's equations where we assume vanishing densities of free charges 
and currents. The dielectric coefficient e{x) and the permeability are periodic and complex, e (a?) = e {x + a), 
fj, {x) = fi(x + a). Here a is any elementary vector of the periodic lattice.^ For given primitive lattice vectors oi 
and 02 the elementary cell 51 C M'^ is defined as = {x E M."^ \ x = aidi + a2a2', < ct2 < 1} x [zsub, Zsup]- A 
time-harmonic ansatz with frequency oj and magnetic field H{x, t) = e^*"*iJ(af) leads to the following equations 
for H{x): 

• The wave equation for the magnetic field: 

W X X H(x) -uj^n(x)H(x) ^0, xen, (1) 

e{x) 

• The divergence condition for the magnetic field: 

V • fi{x)H{x) = 0, xen, (2) 

• Transparent boundary conditions at the boundaries to the substrate (at Zsub) and superstrate (at z^up), 
d^l, where if™ is the incident magnetic field (plane wave in this case), and n is the normal vector on df!,: 

1 



V X (i? - H"') ] xn = DtN{H - ff"), x e dn. (3) 



e(f) 

The DtN operator (Dirichlet-to-Neumann) is in our case realized with the PML method (perfectly matched 
layer). ^ It is a generalized formulation of Sommerfeld's radiation condition and can be realized alterna- 
tively, e.g., by the Pole condition method. 

• Furthermore, the Bloch theorem applies for wave propagation in periodic media. Therefore we aim to find 
Bloch-type field distributions^ solving Equation (1), defined as 

H{x) = e'^'^M(f), u{x) = u{x + a), (4) 

where a is any elementary vector of the periodic lattice and the Bloch wavevector fc G is defined by the 
incoming plane wave if™. 

Similar equations are foimd for the electric field E{x,t) = e~*"*£J(a;). 



For applying the method of finite elements to solve Equations (1) (4), resp. the eorresponding equations 
for the electric field, it is necessary to reformulate these in the weak form. Care has to be taken in the right 
choice of functional spaces which are then discretized by Nedelec's edge elements. This leads to a large sparse 
matrix equation (algebraic problem). For details on the weak formulation, the choice of Bloch-periodic functional 
spaces, the FEM discretization, and our implementation of the PML method we refer to previous works.^'^'^^ 
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Table 1. Geometrical parameters and material parameters for the SRR simulations corresponding to the experimental 
measurements reported by Enkrich et al.^ (See original publication for images with higher resolution.) 

To solve the algebraic problem on a standard personal computer we use either standard linear algebra de- 
composition techniques (LU-factorization) or iterative methods, depending on the problem size. Multi-grid 
algorithms'^ are used within the algorithm for preconditioning. As finite-element ansatz functions, we typically 
choose edge elements of second order. Due to the use of multi-grid algorithms, the computational time and 
the memory requirements grow linearly with the number of unknowns. 

4. SPATIAL DISCRETIZATION AND MATERIAL PARAMETERS 




Figure 3. Visualization of a part of the tetrahedra of the spatial discretization of the SRR geometry shown in Figure 2. 
Dark gray tetrahedra: gold SRR; light gray: air; gray: ITO. Prism elements discretizing the exterior domain are not 
shown. (See original publication for images with higher resolution.) 



Figure 2 shows the geometry of the unit eell of the array of investigated spht-ring resonators. The geometrical 
parameters (corresponding to the experimentally realized SRR, see Fig. 1) are listed in Table 1. This unit cell 
corresponds to the computational domain fl (see Equation 1) of the FEM simulations. 

We discretize the computational domain using a 3D mesh generator. This leads to a coarse unstructured 

tetrahedral mesh containing around 10^ tetrahedra. During the execution of the simulation this mesh is refined to 
a finer mesh (several refinement steps, in each of which every tetrahedron is refined to eight smaller tetrahedra). 
Figure 3 shows elements of a refined mesh. Obviously, by using such unstructured meshes, irregular geometries 
with nearly arbitrary shapes of the scatterers can be resolved ideally. In order to realize transparent boundary 
conditions with the PML method we discretize the exterior space with prism elements above and below the 
sample. The coarse grid in this case contains 1680 prisms which support second-order finite elements. 

In the investigated regime, we assume a dependency of the permittivity of gold on the light frequency given 
by the Drude model: 

e(a;) = l-— , (5) 

with the plasma frequency Up and the collisional frequency lUc- The used parameters as well as the assumed 
permittivities of the other present materials are given in Table 1. 



5. NUMERICAL RESULTS 

The discrete problem corresponding to Equations (1) - (4) with the parameters of Chapter 4 leads to a matrix 

equation with TV = 28 600 unknowns for the coarse grid, resp. A'^ = 133 326imknowns for the one time uniformly 
refined grid. This equation is solved by LU factorization (package 'PARDISO'^^) on a standard 64bit PC {AMD 
Opteron). Typical computation times are 30 sec {N = 28 600), resp. 5min {N = 133 326). 
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Figure 4. Cross sections through the projection onto the x — y— plane of the 3D vectorial solutions at the center of 
the SRR (15 nm above the ITO layer), (a) Electric field, (b) magnetic field. Parameters: A = 1560 nm, incident wave 
polarized in the x — a— plane {E), resp. in the y — a— plane {H), perpendicular incidence. Please note that the magnetic 
field in the inner of the SRR points mainly in the z— direction, which is not seen in this presentation of the data.(See 
original publication for images with higher resolution.) 

Cross-section of the 3D vectorial solutions, E{x) and H{x), for z = const, at the center of the SRR is 
shown in Figure 4. In this case the wavelength of the incident plane wave is A = 1.56 /itm, and the incidence 
is perpendicular (fc = |27r/A| x (0, 0, —1)). The polarization of the electric field of the incident wave is in the 
X — z— plane, i.e., E'o = (1.0, 0, 0) V/m. 

In order to obtain the transmission coefficient for light scattered into the zero*/i diffraction order we perform 
a Fourier transform of the solution at the bottom of the computational domain [z = Zsub)'- 



Aikpc) = / / Ei{x,y,Zsub)e^p{-ikFcx)dxdy (6) 

Clx • dy J-a^/2 J-ay/2 

where kpc is the projection of the wavevector of the zero*/i diffraction order onto the x — y— plane {kpc = 
for perpendicular incidence). In accordance with the convention used in the experiments^'^ we then define the 
transmission as T = ^^7^7 where a corresponds to the transmission of a plane wave through a sample without 
SRR's, lin is the intensity of the incoming wave, and It is the intensity of the transmitted zero-order plane wave 
corresponding to the Fourier coefficients Aj. The reflection, i?, is defined accordingly. 




Wavelength [microns] 



Figure 5. Transmission and reflection of a lightfield incident onto an SRR. Parameters correspond to the experiment 
of Linden et al.,'^ Figure 2 A. Please note the very good quantitative agreement with the experimental results. Trans- 
mission and reflection have been simulated on a coarse grid (x's and solid, resp. dashed, lines) and on a refined grid 
(squares/circles). (See original publication for images with higher resolution.) 

Figure 5 shows the calculated transmission and reflection of a light flcld with perpendicular incidence on a 
SRR. The comparison of the results obtained from simulations on the coarse grid (with iV = 76 014 in this case) 
with the results obtained on a refined grid {N = 374 366) for several wavelengths (squares and circles in Figure 5) 
shows that for the investigated regime simulations with around 10^ unknowns are already well converged and 
the error in T and R can be estimated to be less than few percent. 

The physics of the resonances of the SRR excited by the incident light field (around A = 1.25 /iin and 
A = 3.2 ^m in Figure 5) has been explained in detail in previous works. ^'^ As an example of the strong 
dependence of the resonances on the geometrical parameters of the SRR's we present in Figure 6 two transmission 
spectra for light fields (perpendicular incidence, E— field in the x — z— plane) incident on SRR's with parameters 
according to Table 1, except for the width of the lower bar of the 'U' which is varied by a width of 10 nm. 

By applying light fields which are iiic;idcnt onto the arrays of SRR under oblique angle, it is possible to 
observe several effects: Resonances can be excited now by the magnetic field component of the incident light 
field, resonances are shifted, and new, asymmetric resonances can be observed. Figure 7 shows simulated spectra 
for different angles of incidence. In Fig. 7a, the k vector of the incident plane wave is in the y — 2:— plane and 
encloses angles of a = 0, 30, 45, 60 deg with the z— axis, in Fig. 7b, the k vector of the incident plane wave is in the 
X — z— plane and again encloses angles of a = 0, 30,45, 60 deg with the 2;— axis. The electric field is polarized in 
X— direction (a), resp. y— direction (b) (compare Figure 2 in ref^). All simulations are performed using a refined 
grid with N = 133 326 unknowns. It is again noted that in the case of oblique incidence special care has to be 
taken about the appropriate periodic boundary conditions in x— and y— direction (see Equation 4). In very good 
quantitative agreement with the experimental results^ the main effect of an increased angle a in case (a) seems 





Figure 7. Transmission spectra of light fields incident onto an SRR for different angles of incidence. Polarization of the 
incident field in x— direction (a), resp. j/— direction (b). (See original publication for images with higher resolution.) 



to be an increase in linewidth of the resonance around 1.55 /um. This can intuitively explained by the fact that 
in this case the resonance is excited by a coupling of the electric field to the lower bar of the 'U'-shaped SRR, 

which is not affected by a change of angle a. However, the broadening can probably be attributed to the fact 
that the cross-section of the lower bar has a substantially larger horizontal dimension than vertical dimension. 

The resonance around 1.55 /zm can not be excited by the electric field in case (b). But, as the angle a is 
increased the 2;— component of the magnetic field of the incident wave is increased, too. This component couples 
to the resonance of the SRR. Therefore, around A = 1.55 /im in case (b), and for a = 60deg a weak resonance 
of the SRR can be observed. Again, position as well as strength of the resonance are in quantitative very good 
agreement with the experimental results.^ 

Interestingly, in case (b) the resonance around A = 850 nm is shifted to higher wavelengths for increased 
angle a, which can be attributed to the phase shift between the induced currents in the two upper arms of the 
'U'-shaped SRR. Moreover, another resonance around A = 780 nm appears which is not excited in the case of 
perpendicular incidence of the plane wave. Figure 8 shows cross-sections through the real part of the electric 




Figure 8. Cross sections of the solution in a; — y— direction in the region of the SRR and at the center of the SRR (15 nm 
above the ITO layer). Shown is the projection of the real paxt of the £!— field onto the x — y— plane for different phases 
Lot of the time-harmonic solution. The exciting light field has a wavelength of A = 780 nm and is incident with an angle of 
a = 60 dog (solution corresponds to the low- wavelength resonance in Fig. 7b). (See original publication for images with 
higher resolution.) 

field present in the SRR at different phases of the time-harmonic solution. From these distributions it can be 

seen that the resonance consists in oscillations of the electric field in both upper arms of the 'U', with different 
amplitudes and phases, and an oscillation in the lower bar of the 'U', which leads to an accumulation of charge 
at the outer bottom edges. 



6. CONCLUSION 



In this paper we have investigated rigorous numerical sohitions of the 3D electromagnetic scattering problem of 
plane waves incident onto a periodic array of split-ring resonators in the optical regime. The solutions, obtained 
on standard personal computers, show an excellent agreement with experimental observations. 

Future research directions include the investigation of systems with a macroscopically negative refractive 
index and nonlinear properties of metamaterials. 
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